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ABSTRACT 

The hnear growth factor of density perturbations is generahy beheved to be a powerful observable 
quantity of future large redshift surveys to probe physical properties of dark energy and to distinguish 
among various gravity theories. We investigate systematic effects on determination of the linear growth 
factor / from a measurement of redshift-space distortions. Using a large set of high-resolution A'^-body 
simulations, we identify dark matter halos over a broad mass range. We compute the power spectra 
and correlation functions for the halos and then investigate how well the redshift distortion parameter 
P = f /b can be reconstructed as a function of halo mass both in Fourier and in configuration space, 
where b is the bias parameter. We find that the /3 value thus measured for a fixed halo mass generally 
is a function of scale for k > 0.02 h Mpc"^ in Fourier space or r < 80 h^^ Mpc in configuration 
space, in contrast with the common expectation that /3 approaches a constant described by Kaiser's 
formula on the large scales. The scale dependence depends on the halo mass, being stronger for 
smaller halos. It is complex and cannot be easily explained with the exponential distribution function 
in configuration space or with the Lorentz function in Fourier space of the halo peculiar velocities. 
We demonstrate that the biasing for smaller halos has larger nonlinearity and stochasticity, thus the 
linear bias assumption adopted in Kaiser's derivation become worse for smaller halos. Only for massive 
halos with the bias parameter b > 1.5, the /3 value approaches the constant predicted by the linear 
theory on scales of fc < 0.08 h Mpc~^ or r > 30 Mpc. Luminous red galaxies (LRGs), targeted 
by the Sloan Digital Sky Survey (SDSS) and the SDSS-IH's Baryon Oscillation Spectroscopic Survey 
(BOSS), tend to reside in very massive halos. Our results indicate that if the central LRG sample 
is used for the measurement of redshift-space distortions, fortunately the linear growth factor can be 
measured unbiasedly. On the other hand, emission line galaxies, targeted by some future redshift 
surveys such as the BigBOSS survey, are inhabited in halos of a broad mass range. If one considers 
to use such galaxies, the scale dependence of /? must be taken into account carefully; otherwise one 
might give incorrect constraints on dark energy or modified gravity theories. We also find that the (3 
reconstructed in Fourier space behaves fairly better than that in configuration space when compared 
with the linear theory prediction. 

Subject headings: cosmology: theory — cosmological parameters — galaxies: halos — large-scale 
structure of universe — methods: statistical 



1. INTRODUCTION 

The presence of dark energy, which changes the grav- 
itational assembly history of matter in the universe, ex- 
plains observed acceleration of the cosmic ex pansion well 
within the framework of general relativity (jRies d IT998I : 
iPerlmutter et al.lll999t ISpergel et aT1l2003l . see Komatsu 
et al. 2010 for the latest constraints). There are also 
many attempts to explain the acceleration without dark 
energy by modif ying general relat i vity on cosmologica l 
scales (see, e.g.. iDvali erall [2000l: iCarroll et all [200l . 
Cosmological models in different gravity theories that 
predicts a similar expansion rate H{z), can have the dif- 
ferent cosmic growth rate f{z). The f{z) is often pa- 
rameterized as f{z) — rij^(z) where ^^(z) is the mass 
density parameter at a given redshi ft z and the growth 
index 7 ~ 0.55 in the ACDM model (|Lindedl20Q5l ). Thus 
the precise measurement of the growth rate enables us to 
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investigate the deviation of gravity from the general rela- 
tivity. Recent analysis which focused on such deviations 
using weak gravitational lensing data, cosmic microwave 
background data, and type la supernova data, showed a 
good agreernent w ith the pure ACDM model (e.g., see 
IDaniel et al.ll2010l for the latest work). 

One of the most promising tools to investigate modified 
gravity is redshift-space distortion effects caused by pe- 
culiar velocities in galaxy redshift surveys. In linear the- 
ory an d under the plane-parallel approximation, jKaiseil 
(|1987[ ) derived a formula to relate the observed power 
spectrum of galaxies P^**^ {k, /ik) and the true power spec- 
trum of dark matter Pm\k) through 



p(^\k,^,^) = b\l + p^,lfpj^\k), 



(1) 



where (r) and (s) respectively denote quantities in real 
and redshift space, is the cosine of the angle be- 
tween the line of sight and the wavevector k, (3 is 
the linear redshift distortion parameter related to the 
Towth rate as /3 = f/b, and b is the bias parameter 
Kaiseii [l984[ ). Thus the measurement of the redshift- 
space distortions allows one to directly probe deviations 
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from general relativity, although the determination of 
the biasing is another important issue. The Kaiser's 
formula (equation (H}) is modified on small scales be- 
cause the nonlinear random velocities smear the clus- 
tering along t he line of sig ht known as the 'finger- 
of-god' effect (jJacksonI Il972[ ). However such a non- 
linear model still relies on Ka iser's formula on large 
scales (jPeacock fc Doddi |1994[) . For the im portance 
of nonlinearity on s u ch scales, see iScoccimar ro (2004) 



Taruva et al.' (200% iDesjacques fc ShethI (|201Q ). and 



Jennings et al. (201C 

Constraints on have been repo rt ed in vari- 
ous surveys (e.g.. Peacock et al.l |200H iZehayi et al.l 
200l Hawkin sgball l2003t iTegmark et al.l 12004 120061: 
Ross et al. .20071: i Guzzo et al.' '20 081: Ida Ang ela et ahl 
2008; Cabre_&£aztaiiaga 2009). O kumura et al. (2008) 
also showed using a luminous red galaxy sample from the 
Sloan Digital Sky Survey (SDSS) that simultaneously 
analyzing redshift-space distortions and anisotropy of 
the baryon acoustic scales allows one to give a strong con- 
straint on dark e nergy equation-of-sta t e, as was theoret- 
ically predicted fHu fc HaimanI 120031 : ISeo fc Eisensteinl 
120031 : Matsubara 2004 ) and th i s fact w as explicitly 
emphasized by fAmendola et al.l (|2005l ). IGuzzo et al.l 
(|2C)08i) considered constraints on / to test the deviation 
from general relativity us ing the observations at dif- 
ferent redshifts (see also, 'Di Porto fc Amendola' '2008': 
Nesscris & Pcrivolaroooulos 2008; Yamamoto ct al. 
2008t ICabre fc Gazt aiiaga 2009). We note that all these 



previous studies have used linear theory prediction of 
redshift-space distortions to compare with their measure- 
ments on scales presumabl y large enou g h for the linear 
theory to be valid. Na kamura et al.l ()2009[ ) adopted 
an alternative approach and constrained the growth 
factor by m easuring the dampi ng of the baryon acoustic 
oscillations. iReves et al.l (|2010D gave a strong constraint 
on modified gravity theory and confir med general 
relativ ity using the method proposed by iZhang et al.l 
((200l which can eliminate the uncertainty of the galaxy 
biasing by combining weak gravitational lensing, galaxy 
clustering, and redshift-space distortio ns (for similar 
theoretical attempts, see e.g.. iPercival fc Whitd 120091 : 
IMcDonald fc Seliakll200l ISong fc Kavdl2010n . 

There are many ongoing and upcoming large 
galaxy surveys, such as the SDSS-IIF s Baryon Os- 
cillatio n Spectroscopic Survey (BOSS; iSchlegel et al.l 
2009al). the F iber Muhiobject Spectrograph (FMOS; 



Sumivoshi et a l. 2009), t he Hobby-Eberl y Dark Energy 



Experiment (HETDEX: iHiU et al.l Hool the WiggleZ 
(IGlazebrook et al.l 120071) . the BigBOSS (ISchlegel et al.l 
l2009bD . and" so on. These observations are expected to 
enable us to distinguish among gravity theories with high 
precision through measurement of redshift-space distor- 
tions as well as that of baryon acoustic oscillations. How- 
ever, it is not clear if the accuracy of predicting redshift- 
space distortions is better than or comparable to the pre- 
cision required from future surveys. In addition, we do 
not know how large the deviation from true cosmology is 
if any. Precision of the constraint may depend on galaxy 
types, such as luminosity and host halo mass. There 
were many attempts to investigate the validity to use the 
redshift-space distortions to extract the cosmological in- 
formation (e.g.,'Hatton fc Colen998',l^;'Beiiind_eLalJ 
[2001; .Tinker et al...200& ). .Tinker et al. (2006. ) found that 



/3 can be estimated accurately using linear theory if the 
finger-of-god effect is removed perfectly. 

In this paper, we present a detailed study on this aspect 
using a large set of iV-body simulations. We measure the 
power spectra and correlation functions of dark matter 
halos. Using them, we estimate the redshift distortion 
parameter /? from the monopole-to-real-space ratio and 
the quadrupole-to-monopole ratio, both of which are re- 
lated to /3 in linear theory. Then we examine whether 
l3 measured in these ways can give true cosmological in- 
formation. We also investigate the dependence of the 
precision of reconstructed /3 on halo mass. Particularly 
we will clearly show that the /3 value obtained from the 
small-halo clustering does not approach a constant even 
on large scales as linear theory predicts. In addition such 
small halos are shown to be more stochastic tracers of 
the underlying density field than massive halos. We also 
discuss in detail on which scale and with which method 
one can get the correct (3 or f from the redshift-space 
distortions. 

The paper is organized as follows. In Section [2l we 
describe the A^-body simulations and the halo occupa- 
tion distribution model used to populate them with mock 
galaxies. The basic two-point statistics used in our anal- 
ysis are also presented. In Section [3] we briefly review lin- 
ear theory of redshift distortions and how to determine 
the redshift distortion parameter /? from the power spec- 
trum and the correlation function. Section [J] is devoted 
to the analysis of redshift distortion effects in simulations 
to determine /3 and the growth rate /. Our conclusions 
are given in Section [5] 

2. Af-BODY SIMULATIONS 

2.1. Dark Matter Halo and Galaxy Catalogs 

We use a large s et of A^-body si mulations, which is an 
updated version of lJing et al.l (j2007il , to create dark mat- 
ter halo distribution. We assume a spatially flat ACDM 
model with the mass density parameter flm — 0.268, 
the baryon density parameter fib — 0.045, and the Hub- 
ble constant h — 0.71. Initial conditions are gener- 
ated using the matter transf er function by CMBfast code 
(jSeliak fc Zaldarriagalll996t ) and the density fluctuation 
amplitude is set to be as = 0.85. We employ 1024^^ 
particles in 15 cubic boxes of side 600 Mpc and 4 
of side 1200 Mpc, respectively abbreviated to L600 
and LI 200. We mainly show results obtained from the 
L1200 boxes, while the L600 boxes are used in order to 
analyze dark matter halos with small mass and to see if 
the L1200 samples have good enough resolution for the 
smallest halos. Simulation parameters are summarized 
in Table [TJ See Jing et al. (2007) for details of the simu- 
lations. Dark matter halos are identified at the redshift 
•^out using the friends-of-friends algorithm with a linking 
length equal to 0.2 times the mean particle separation. 
All unbound particles in the EOF halos are further dis- 
carded. As shown by I Jing et al.l (|2007D . it is necessary 
to eliminate these unbound particles in order to have a 
correct measurement of clustering for small halos of a 
few tens particles. We use all the halos with more than 
12 particles. Identification of small halos is subtler than 
that of massive halos because of the limited number of 
dark matter particles which constitute small halos. As 
will be shown with L600 and L1200 simulations in Section 
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TABLE 1 

Simulation parameters 



boxsize 


particles 


realizations 




^out 


600 


10243 


15 


1.5 X 10"' 


0.295 


1200 


10243 


4 


1.2 X 10^1 


0.274 



Note. — nip in column 4 is the particle mass. 



TABLE 2 

Properties of simulated halos and galaxies 



box 


M(h-^MQ) 


rip 




bik) 


b{r) 


L600 


2.2 X 10" 


13 < np < 18 


1.3 X 10"^ 


0.69 


0.70 




1.7 X 10^2 


92 < rip < 136 


1.9 X 10^ 


0.88 


0.89 


L1200 


1.8 X 10^2 


12 < rip < 17 


1.7 X 10'^ 


0.89 


0.88 




1.4 X 10^3 


92 < rip < 136 


2.2 X 10^ 


1.30 


1.30 




1.0 X lO^-* 


692 < rip < 1037 


2.2 X 10* 


2.28 


2.31 




LRG 


12 < np < 25000 


1.4 X 10^ 


1.90 


1.94 















Note. — The halo mass M in column 2 shows the central 
values of each mass bin. Column 3 shows the ranges of the number 
of particles. Afjialo is the total number of halos. b(k) and b{r) in 
columns 4 and 5 are the best fit bias para mete rs in Fourier and 
configuration space, respectively (see Section l2.2l l. 



I2.2l and Section|4l however, the clustering of halos can be 
well measured to this limit. We choose Zout ~ 0.28 be- 
cause the luminous red galaxies (LRGs) of the SDSS are 
at this redshift, but almost all of our conclusions should 
not rely on our choice of this particular redshift. 

We consider as a m ock g alaxy catalog the LRG sam- 
ple jEisenstein et al.l 120011 ) in the SDSS ()York et al.l 
120001 ). In order to populate the center of the halos 
with LRGs, we rely on the fra mework of the halo oc- 
cupation distribution (HOD, e.g.,'JinE et al."1998VSeliay 
2000: Scoccmiarro et al. 2001; Bcrhnd & Weinberg 2005 



Yang et al.H2003l: IZheng et al.ll2005D , which describes the 



relationship between the galaxy and dark matter density 
fields. Galaxies are assigned to the hal os using the best 
fit HOD parameter s for LRGs found bv lSeo et all (pOOl ) 
(see also IZheng eT al. 2005, 200 ^. This meth od was ap- 
plied in our previous work (jOkumura et al.l [2009) and 
the good agreement with the observation in clustering 
has been confirmed. LRGs are found to reside in mas- 
sive halos of typical mass ~ 2 x 10" - lO^/i-^M©. The 
fraction of satellite LRGs is 6.3% and only central LRGs 
are used for our analysis below. The peculiar velocity of 
their halos is assigned to central LRGs. Table [5] lists the 
detail of the representative halo and LRG catalogs. 

2.2. Two-point Statistics 

We plot the two-point statistics for dark matter halos, 
central LRGs, and dark matter particles measured from 
our simulation catalogs in Figure [T] When dark matter 
halos are analyzed, they are divided into narrow mass 
bins as Mi - AMi < Mi < Mi + AMi for ith bin where 
AMi = 0.2Mi. Corresponding ranges of the number of 
particles for each halo catalog are listed in Table [2] Re- 
sults of LRGs and dark matter are shown only for the 
L1200 samples while those of high and low mass halos 
are shown for the L1200 and L600 samples, respectively. 
Results are averaged over all the realizations. In redshift 
space, positions of objects are displaced as a result of the 



peculiar velocity along the line of sight. We regard each 
direction along the three axes of simulation box as the 
line of sight and the statistics are averaged over three 
projections of all realizations for a total of 45 samples 
for the L600 simulation and 12 for the L1200 simulation. 
The error bars shown in Figure [T] are the standard error 
of the mean. 

Figures [1] (a) and (b) show the power spectra in real 
space, P('')(fc), and in redshift space, p('')(fc), respec- 
tively. In both figures clear halo-mass dependence of the 
clustering amplitude can be found. The gray and green 
lines show the halo power spectra with the same halo 
mass in the L600 and L1200 samples, respectively. It can 
be easily seen that the agreement of the power spectra be- 
tween the two boxes is very good for k > 0.03 h Mpc~^, 
indicating that resolution of a halo with > 12 particles is 
enough for the clustering analysis here. The discrepancy 
between the two lines at the smaller k < 0.03 h Mpc~^ 
is owing to the large cosmic variance effect in the L600 
sample. Figures [1] (c) and (d) show the correlation func- 
tions in real space, ^^''•'(r-), and in redshift space, ^'■-^^r), 
respectively. The suppression of the correlation func- 
tions for halos and central LRGs is caused by the finite 
size of the halos and it is alleviated to some e xtent in 
redshift space due to their random velocities (jJacksoni 
119721) . The baryonic acoustic features, which were clearly 
detected in our simulation samples as a single peak at 
~ 100 Mpc, appear as slight bumps in Figures [1] (c) 
and (d) because we use the binning much broader than 
the width of the peak. 

The bias parameter can be computed both in Fourier 
space, b{k), and in configuration space, b{r), through 



bik) 



>W(fc)' 

, Pm "* (k) ^ 



1/2 



1/2 



(2) 



Figure [2] shows the b{k) and b{r) for dark matter ha- 
los and LRGs. The results are averaged over realiza- 
tions and the error bars are the error of the mean. In 
both configuration and Fourier space, the bias parame- 
ters for halos and LRGs are found to be almost constant 
on sufficiently large scales. We assume the bias to be 
constant and search for the best fit value for each sam- 
ple by computing the statistics. We compute for 
16 < r < 79 h^^ Mpc in configuration space and for 
0.018 < k < 0.10 h Mpc~^ in Fourier space for results 
from the L1200 simulations. On the other hand we com- 
pute it for 16 < r < 63 Mpc in configuration space 
for the L600 samples because the large-scale data is not 
very reliable owing to the cosmic variance while we still 
use the same range in Fourier space. Figure [3] shows 
the halo biasing as a function of the halo mass, b{M). 
The error bars show the 95% confidence interval. The 
results between L600 and L1200 are largely overlapped 
with each other and we confirm that the systematic er- 
ror caused by the different box sizes is negligibly small. 
In addition, the bias parameters obtained in Fourier and 
configuration space are consistent. The best fit values 
obtained here are used for the theoretical prediction of P 
through Kaiser's formula in Section |4l 
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Fig. 1. — Two-point statistics for dark matter halos, LRGs, and dark matter particles, (a) Real-space power spectra P^^\k). (b) 
Redshift-space power spectra P(^'(A;). (c) Real-space correlation functions ^''^^(r). (d) Redshift-space correlation functions ^^^\r). The 
values quoted in the figure are the halo mass at the center of mass bin in units of the solar mass [/i~^Mq]. Error bars are the standard 



error of the mean. 

There are at least two ways to deterniine the red- 
shift distortion parameter /3. They have been well de- 
velop ed both in Fourier space (iKaiseii 119871: iCole et al.l 
I1994D and in configuration space (|Haniiltoi]| |1992') under 
the plane -parallel approxi mation and summarized in a 
review bv lHamiltonI (|1998D . who also collected the obser- 
vational constraints then available on /? in various sur- 
ve ys. We follow the si milar notation with that adopted 
bv [Tinker et"all pQOl ). 

3.1. Fourier Space 

For plane-parallel redshift-space distortions, t he red- 
shift space power spectrum can be written as (jKaiseri 
[19871) 

(3) 

where Li are Legendre Polynomials. The multipoles of 
the redshift-space power spectrum are expressed as 



real-space power spectrum p(°/'') and the quadrupole-to- 
monopole ratio P^^/^^ (iCole et al.lll994| ): 



P 



Pojk) 

p(^Hk) 

P2{k) _ 
Po{k) 



i + -/3 + -r 



1/3 + 1/3' 



1 



(5) 
(6) 



Piik) = 



21 + 1 



+1 



F(^)(fc,Aik)i/(A'k)dMk- (4) 



The last equality in the two equations holds only on large 
scales where linear theory can be applied. 

3.2. Configuration Space 

The redshift-space correlation functions can be ex- 
pressed similarly to the power spectra under the plane- 
parallel approximation as 

i^'Hrp,r^) ^ Ur)Lo{^l)+^2ir)L2{^i) + Ur)L4{^i), (7) 

where rp and r^r are the separations perpendicular and 
parallel to the line of sight and fi is the cosine of the 
angle between the separation vector and the line of sight 
/i — cos 9 = r^r/r. The multipoles of the redshift-space 
correlation function are expressed as 



We can derive two useful combinations of them which 
are directly related to /3, the ratio of the monopole to the 



6(0 



21+1 



^^'\rp,r^)Li{ij)dii. 



(8) 



Determination of the Growth Factor from Redshift Distortions 



5 



2 - 



1 - 



' 1 1 1 1 1 1 1 1— 

4. i (a) 

^fA.....j«i.....A,....A,....^....^. 




1 - 



— I — I — I — r- 



(b) 



i6t 



:©■-€>-©. 



■€>■© 



»-o-b-b-bH-b-b-»- 



_l I I I I 



— I 1 1 1 1 — I — 1 — r 



I I l_ 



_l U I I L_ 



0.01 



0.1 



10 



100 



k[h Mpc '] 



r[h"^Mpc] 



Fig. 2. — Halo bias (a) in Fourier space and (b) in configuration space. The color of each line corresponds to the one with the same 
color in Figure [T] The bias for the central LRGs is also plotted for comparison. Error bars are the standard error of the mean. The region 
enclosed by the two vertical lines shows the scales where we assume the scale-independent bias. 



2 - 



1 - 



— 1 — I I I I 1 1 



1 1 — I I I I 1 1 1 

P(k) LI 200 

L600 

^(r) LI 200 

L600 



—I — I I I I I I 




10 



11 



^12 



^13 



10'^ 10"" 10 

Halo Mass M [h'^M®] 



14 



Fig. 3. — Halo bias as a function of the halo mass. Shown is the 
best fit values for constant fits of b between 15 < r < 80 Mpc 
(L1200) and 15 < r < 60 h'^ Mpc (L600) for C(r) and 0.02 < fe < 
0.1 h Mpc"^ for P(fc). Results from the L1200 samples are shown 
for 1.4 X 10^2 < M < 1.9 X lO^^M© and those from the L600 
samples for 1.8 X lO^^M < 1.6 X IO^^Mq. The error bars show the 
95% confidence levels. 

In hnear theory, the ratio of the monopole to the 
real-space correlation function and the quadrupole-to- 
monopole ratio are relate d to the redshift distortion pa- 
rameter 13 on large scales ()Hamiltonlll992[ ): 



Coir) 



2 1 , 
= l + -/?+-/32, 



(2/0) ^^W. 



I/?- 



Coir) -Coir) 



1 



13 



(9) 
(10) 



where ^o('") — i-^/r^) Jq Coir')r'^dr' . When one wants 
to constrain the pairwise velocity dispersion of galax- 
ies which becomes dominant on small scales, the real 
space correlation function is convolved with the distri- 
bution function of pairwise velocities to g ive the redshift 
space correlation function (jPeeblesI 11980), which is not 
the purpose of this paper (s ee, e.g .. [Peacock et al."2001'; 
iZehavi et all [2002; .Hawkins et al.l [2003; .Jing k Borner, 



[lMlGuzzo"eran[2rMrCabre fc Gaz tafiaga"2009y We 
will briefly discuss the effect of the pairwise velocities on 
f3 reconstruction at Section 14.21 

4. RESULTS AND DISCUSSION 
4.1. /3 Reconstruction 

In Figure [4| we show the resulting /3 values of dark 
matter halos, LRGs, and dark matter reconstructed by 
the methods described in Section [3] In each panel the 
horizontal lines show the large-scale values predic ted by 
general relativity, /3 = n°J^iz)/b (jLindeii |2005[) . For 
the bias parameters in Fourier and configuration space 
we use the best fit values obtained in Figure [3] The (3 
value of dark matter is simply equal to the growth rate 
/ because 6 = 1. We can see the agreement of the f3 val- 
ues obtained from the L600 and L1200 samples with the 
same halo mass, thus the different number of particles, 
indicating that resolution of a halo with 12 particles is 
accurate enough for /3 reconstruction. The discrepancy 
between the two on large scales is again owing to the 
cosmic variance in the L600 sample. 

Figure [4| (a) shows /3 as a function of k measured from 
the ratio of the monopole to the real-space power spec- 
trum 

p(o/r) ^p^/p(r)^ Small-scale values obtained from 
the dark matter particles and small-mass halos are sup- 
pressed by the random peculiar velocities, while those 
from the LRGs and massive halos go up on such scales 
as a result of the finite size of halos, as described in 
Section 12.21 One can see that, for the most massive 
halos (M - IQI^/i-IMq) and LRGs, the ratio P(°/'^) 
reconstructs the /3 values predicted by linear theory at 
k < 0.08 h Mpc"^ This means that linear theory is 
accurate enough to predict and constrain /3 using such 
massive halos. On the other hand, the smaller halos we 
focus on, the more prominent scale dependence of /3 we 
find even on large scales. This behavior of /3 from the 
smaller mass halos is consistent with that obtained by 
[Tinker et al.l ([5006). They used mock galaxies assigned 
to the halos by t he HOD parameter s applied to the SDSS 
MAIN galaxies ([Zehavi et al.[ [20051 ) which preferentially 
res ide in halos w ith small mass. Thus the behavior seen 
bv [Tinker et al.i (,2006.) is found to be caused by the con- 
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Fig. 4. — Redshift distortion parameter /3 reconstructed from; (a) the monopole-to-real-space ratio of the power spectra; (b) the monopole- 
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color and line type, where the best fit parameter for the biasing is used for the prediction. Error bars are the standard error of the mean. 
The diamonds and the open circles have been offset in the horizontally positive direction for clarity while the open squares and the triangles 
in the horizontally negative direction. 



tribution from smaller halos. Note that here we used 
the real-space power spectrum measured from the sim- 
ulations. It is not a direct observable and translation 
of the power spectrum from redshift space to real space 
usually draws the additional error. 

Figure S] (b) is the same as Figure S] (a), but the /3 
values are measured as a function of separation r com- 
puted from the ratio of the redshift-space and real-space 
correlation functions ^("Z'') — ^a/^^'^\ The slight differ- 
ence of the linear theory prediction between Fourier and 
configuration space is due to the diff'erence of the best fit 
parameters for the biasing seen in Figure [S The scale 
dependence of the reconstructed (3 found by using ^("Z'') 
is more prominent than that by using P^^/^K Even the 
result obtained from the LRGs is monotonically increas- 
ing, intersects with the prediction from linear theory, and 
does not draw closer to a constant on all scales probed. A 
similar beh avior has also been found for dark matter by 
iCabre fc~Ga ztafiaga (2009) up to 40 Mpc (see also 
the red line in Figure 4b), but the tendency is much 
more significant for dark matter halos, even for those 
with b ^ 1. Constraints on /3 are usually given under 
the assumption of one constant parameter over a scale 
range for which the statistics is computed. However, 
according to Figure |4] (b), it could be a coincidence that 



one gets the true value of /3 as a best fit parameter. Thus 
one should be cautious when the deviation from general 
relativity is investigated through the measurement of /3 
from the ratio S.o/S,''^^- 

In measuring /3 from the quadrupole-to-monopole ratio 
in Fourier space P^^/*') = P2/P0, one needs to measure 
P^^\k,ii) in finite bins, usually constant separations in 
fi, and numerically integrate it along fi direction. Hence 
the finite bin size may cause a systematic error in mea- 
surement of /3. Using linear theory, we test the accu- 
racy of the integration between constant fi and constant 
9 = cos~^ fi binnings. We found constant fi binning un- 
derestimates (3 by 2.5% while constant 6 binning overes- 
timates by 1.3% for 10 bins. We thus adopt the constant 
9 binning and take the number of bin to be 10 between 
< 9 < 90°. Figured] (c) shows /3 measured from the 
quadrupole-to-monopole moments p(-^/'^) . This quantity 
can be directly measured in observation. We put artifi- 
cial large-scale cuts in (3 values measured from the L600 
samples because they have limited number of modes and 
thus are strongly affected by the cosmic variance. We 
can see this by the difference between the green and gray 
lines because they have similar halo mass. On the other 
hand, we did not use such strong scale cuts in Figure |4] 
(a) because the effect of the cosmic variance can be clim- 
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Fig. 5. — Quadrupole-to-monopolc ratio in Fourier space for dark 
matter particles (red) and for massive (magenta) and small (light 
blue) halos. The horizontal lines are the linear theory prediction. 
The upper and lower dashed curves in red are the non-linear predic- 
tion with the pairwise velocity dispersion cr^ = 500 and 600km/s, 
respectively, while those for halos are with cry = 400 and 500km/s. 

inated t o some extent by taking t he ratio of two power 
spectra (jMcDonald fc Seliak|[2009l) . In Fourier space, the 
behavior of 13 measured from p(2/°) is almost the same 
as that from p(o/'') except for the magnitude of the error 
bars owing to the cosmic variance. 

Finally, we show /3 values measured from the 
quadrupole-to-monopole ratio in configuration space 
^(2/0) _ ^2/^0 in Figure 2] (d). In measuring the 
quadrupole moment by equation ([S]) we adopt the bin- 
ning on a polar grid of logarithmic spacing in r and linear 
spacing in angle, then numerically integra te the corre- 
lation function at each r, as was done by iTinker et al.l 
([2OO6). On small scales, the behavior of /3 thus de- 
termined is complicated. For small halos, the results 
give lower /? values than those from linear theory pre- 
diction, which may be caused by the random peculiar 
velo cities and thus might be correctable to s ome extent 
(see lHawkins et al.ll2003t iTinker et all l2006( ) . However, 
note th at peculiar v elocities predicted from a simple halo 
model (jYang et al.i 2003) have very different luminosity 
dependence with tho se from observations (i Ji ng fc Sorneii 
l2004HLi et alll2006l see also Slosar et al. 2006). For the 
massive halos, the /3 values on small scales are larger than 
the linear theory prediction, which might indicate that 
these halos are approaching with each other. Interest- 
ingly, for the most massive halos (1.0 x 10^*h~^ Mq) and 
the LRG sample, we can simply use linear theory and use 
the data points on scales larger than ^ 25 Mpc in 
order to constrain the growth rate. For galaxies within 
the halos with the mass ^ 1.4 x IO^-^H^^Mq, the /? value 
is coincidently a constant, but it is lower than the linear 
theory prediction. If we use a population of such galax- 
ies, we would underestimate the growth factor. Lastly 
for galaxies within less massive halos, the reconstructed 
P becomes a scale-dependent function, and one has to be 
extremely careful in extracting the linear growth factor 
from a measurement of redshift-space distortions of such 
galaxies. 

4.2. Pairwise Velocity Dispersion 



In order to see if random peculiar velocities can cause 
the deviation from the linear theory prediction, we con- 
sider a simple Exponential model for the pairwise veloc- 
ity dispersion (PVD) in configuration space which results 
in a Lorentz damping factor in Fourier space, 

G{k,fXk,c7,) = il + eiilal/2)-\ (11) 

The power spectrum of this dispersion model is expressed 
as Equation ([T]) multiplied by G(A:, /^k, cr„). Although 
more accurate models have been developed by many au- 
thors, this simple model is useful enough for our purpose. 
In Figure m we show the predictions for the quadrupole- 
to-monopole ratio in Fourier space p(^/°) from the dis- 
persion model. In order to avoid making the figure un- 
clear, we show the results for only dark matter particles, 
the most massive halos and the smallest halos. The over- 
all shape of p(-^/") for dark matter is well explained by 
the dispersion model with ay ~ 600 km/s, which has 
already been foun d using more accurate models (e.g., 
iTaruva et al.ll2010f ). However, the results for dark mat- 
ter halos are much more complicated. Here we adopt 
ay ^ 450 km/s for the halos according to (jHamana et al.l 
|2003[ ). For small halos, not only the deviation of the 
measurement from the linear theory prediction but also 
its scale dependence cannot be corrected by the model. 
On the other hand, the behavior of /3 reconstructed for 
massive halos is opposite to the dispersion model. These 
complex results are somewhat expected because the halo 
peculi ar velocities change very mildly with their mass 
(jHamana et al.ll20C)3[ l. Thus the difference of the scale 
dependences of /3 among small and large halos cannot be 
simultaneously explained with such analytical models for 
the PVD. We note th at the model of PVD was adopted 
bv lGuzzo et al.l (|2008( ) to correct for the nonlinear effect 
when they measure the /3 at redshift 0.8. 

4.3. Nonlinear Stochastic Biasing 

In the derivation of Kaiser's effect, the linear bias rela- 
tion between objects considered (galaxies or halos) and 
dark matter was assumed. However, there could be con- 
siderable stochasticity and nonlinearity in the bias rela- 
tion between dark ha los (or galaxies) and dark matter 
(|Dekel fc Lahavl I1999D. We follow the formalism pro- 
posed by iTaruva fc Sutol (I200Q) and applied to simula- 
tion data bv lYoshikawa et al.l (I2001D. We briefly summa- 
rize some parameters defined by ITaruva fc Sutol (j2000D 
to quantify the nonlinearity and stochasticity of the halo 
bias. First, the density fields of dark matter and dark ha- 
los are evaluated as 6„i{x, Rs) and (5(x, Rs), respectively, 
smoothed over the top-hat window radius Rs- The bias 
parameter and the correlation coefficient are respectively 
defined by 

Note that the bias parameters defined in Section [2?2] are 
from the two-point statics, while the bias 6var defined 
here is from one-point statistics. In order to quantify the 
nonlinear and stochastic nature of the biasing separately, 
two more useful parameters are introduced. For this pur- 
pose, let us define the conditional mean of 6 for a given 
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where P{S\dm) is the conditional probability distribution 
function. Then the nonlinearity of the biasing is quanti- 
fied by 



{ss„,y 



which vanished only when the biasing is linear. Similarly, 
the stochasticity of the biasing is characterized by 



Fig. 6. — Joint probability distributions of overdensity fields for dark halos with dark matter overdensity in real space {left panels) and 
in redshift space (right panels) smoothed over Rs = 20 Mpc and Rs = 30 Mpc. Upper, middle, and lower panels in each set 
show the results for halos with M = 1.9 X 10^^, M = 1.4 X lO'^^, and M > 1.9 X 10^^ , respectively. The values in the color bars show the 
logarithm of the probability to base 10. The solid lines are the conditional mean S{5m)- The dashed lines show the linear bias relation 

5(l5m) = bvavSm- 

S„i , In Figure [U we show the joint distribution of 6 with 

Sm in real space (left) and in redshift space (right). From 
the top to bottom, the results obtained from the smallest 
mass bin of the L600 samples {M = 1.9 x lO^/i^Ufo), 
the smallest mass bin of the L1200 samples {M — 1.4 x 
IO^'^H'^Mq), and the largest mass bin of the LI 200 sam- 
ples {M > 2.9 X 10^^/i"^Mq) are plotted. In each panel 
we also plot the conditional mean relation 6{6m,) as the 
solid line and the linear bias relation S = &var<5m as the 
dashed line, both of which are obtained from our simu- 
lations. Here we focus on halo mass dependence of the 
nonlinear stochastic biasing. The deviation of 5 from the 
linear bias is caused by the nonli near stochastic bias as 
well as the halo exclusion effect (jYoshikawa et al.ll200lL 
see also Smith et al. 2007). The halo exclusion effect 
is alleviated in redshift space due to the random veloc- 
ities of halos, as we have seen in Section 12.21 Despite 
the fact that this exclusion effect is more significant for 
larger thus more massive halos, the deviation from the 
linear relation for such halos is much smaller. This indi- 
cates that more massive halos have the smaller nonlin- 
earity and stochastic ity, and the latter was also found by 
iHamaus et al.l ()2010[ ) using a complementary statistics in 
real space. 

In order to see these effects more quantitatively, we 
show Tcorr, Eni, and egcatt as functions of halo mass in 
Figure [7] from the top to bottom. The results for the 
parameter Cni show that the nonlinearity of the halo bi- 
asing is smallest for the most massive halos. It increases 
for smaller halos and gets close to a constant. Simi- 



^scatt 



(15) 



which vanishes for the deterministic bias where S = 

Following the same procedure of Yoshikaw a et al.l 
(|2001i) . we compute the parameters defined above for 
our halo catalogs in real and redshift space. In order 
to minimize the Poisson noise effect and fairly compare 
the results of different mass, only in this subsection we 
keep each subsample of a given halo mass to have the 
same number density, 1.16 x 10~'*(/i^^Mpc)~^. This den- 
sity corresponds to A'haio ~ 2.5 x 10** for L600 samples 
and A'^haio ~ 2.0 x 10^ for L1200 samples. We adopt the 
smoothing scales Rg = 20 and 30 Mpc. Many pairs 
of the values [5(x), Jm(x)] are obtained for randomly se- 
lected points x in the simulation box. 
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larly, the stochasticity parameter egcatt has the minimum 
value for the most massive halos. At a whole mass range 
probed, however, the stochastic bias monotonically in- 
creases toward the lowest mass. Finally the stochasticity 
of the smallest halos becomes five times larger than that 
of the most massive halos. Thus both the nonlinearity 
and stochasticity of the halo biasing, particularly the lat- 
ter, is likely one of the cause for the systematic deviation 
of /3 values from the theoretical prediction. 

4.4. Halo Mass Dependence of Growth Rate Constraints 

Because we also calculated the power spectra and cor- 
relation functions for dark matter particles in the same 
samples as those for dark matter halos and LRGs, we 
can directly reconstruct the growth rate through 



/(r) = fe(r)/3(r) = 



/(fc)=6(fc)/3(fc): 



pW(fc) 



1/2 



/3(r) 



1/2 



(16) 



(17) 



The difference between the use of /? and / is just whether 
the bias b is used as prediction or measurement. When 
/ is used, however, we can take into account the slight 
scale dependence of the bias seen in Figure [2] Fig- 
ures [8] show the growth rate measured from the four 
methods described above. The horizontal line shows 
the input ACD M model pred icted from general relativ- 
ity / = n°n^^iz) (|LindeHl2005l ). We can also see that the 
linear redshift distortions for the LRGs reconstruct the 
true value of / well, except for that from ^(o/'"). 

Here let us discuss which method can be used to obtain 
the growth factor / better. Figure M shows the compari- 
son among the 4 methods for the measurement of / from 
the LRG clustering. As we have seen above, we find 
strong scale dependence of the growth factor obtained 



from ^(o/''). On the other hand, both of the methods 
in Fourier space, p(°/'') and pf^/'^), and the estimator 
^(2/0) configuration space, can give good estimation of 
(3 in linear theory on scales A = ir/k > 30 Mpc or 
r > 30 h-^ Mpc. 

5. CONCLUSION 

We have investigated how accurately the redshift-space 
distortions can be used to measure the linear growth fac- 
tor /. The growth factor is a powerful observable tar- 
geted by future large redshift surveys to probe dark en- 
ergy and to distinguish among different gravity theories. 
For this purpose, we constructed a large set of A'^-body 
simulations, dividing each dark matter halo catalog into 
the subsamples with narrow mass ranges. As an example 
of a galaxy sample, mock SDSS LRG catalogs were con- 
structed by applying the HOD modeling to the simulated 
halos. Then we have measured the two-point statistics, 
power spectra and correlation functions, for dark matter 
halos and LRGs. The dark matter halos were analyzed 
as a function of halo mass in order to see dependence of 
the /3 measurement on the halo mass. 

We have determined the /3 values as a function of halo 
mass and scale using four methods. First, we found that 
/3 reconstructed from the ratio of the monopole to the 
real-space power spectra p(°/'") = pi^) I pi^) (equation 
([5])) asymptotically approaches the true value. In partic- 
ular for the massive halos and LRGs, the prediction from 
linear theory known as Kaiser's formula is applicable to 
give a correct constraint on the growth rate. However, 
for less massive halos, the ratio approaches the true value 
only at a very large scale k < 0.02 h Mpc^^. Second, 13 
reconstructed from the ratio of the monopole to the real- 
space correlation function ^(o/'') = ^{^)/^(^) (equation 
^) approaches neither the true value nor a constant 
even on large scales. This statement is valid especially 
for small halos with the bias parameter 6 < 1. Because 
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Fig. 8. — Same as Figure|4]but for growth rate /, reconstructed from; (a) the monopole-to-real-spacc ratio of the power spectra; (b) the 
monopole-to-real-space ratio of the correlation functions; (c) the quadrupole-to-monopole ratio of the power spectrum; (d) the quadrupole- 
to-monopole ratio of the correlation functions. The horizontal line shows the linear theory prediction, Q!^^{z). Error bars are the standard 
error of the mean. The open squares and the open circles have been offset in the horizontally negative direction for clarity while the 
triangles in the horizontally positive direction. 

as but with lar ger error bars as expected. Finally, 

we found that when the quadrupole-to-monopole ratio 
in configuration space ^(2/°) = (equation ([10])) is 

used, a similar conclusion is reached to that of P'^/'^) 
when r = A = 7r/fcis adopted. 

For small halos with b < 1.3, the reconstructed /3 val- 
ues do not approach a constant in most of measurable 
regions, particularly those from ^('^Z'") in the configura- 
tion space. No method can provide a reliable estimator 
for the determination of the growth factor from the clus- 
tering of such small halos on the large range of scales 
probed. Using the halo catalogs with different box sizes, 
we confirmed that such a behavior is not caused by the 
resolution effect of small dark matter halos. While the 
scale dependence changes with the halo mass, the pe- 
culiar velocity of halos do es not change much with the 
mass (jHamana et al.l[2003| ). Using the simple dispersion 
model, we demonstrated that the different scale depen- 
dence of /3 among small and large halos cannot be simul- 
taneously explained. Also there are two types of velocity 
biases which affect the redshift distortion; the dynami- 
cal bias caused by dynamical friction and the spatial bias 
caused by the difference between the distribution of halos 
and that of dark matter. There is no dynamical velocity 
bias because the halo velocities are determined from the 
mean velocities of dark matter within the halos in our 




rorn/k [h'^Mpc] 



Fig. 9. — Comparison of the growth rate measurement from LRGs 
among the different methods. The horizontal axis shows separation 
r for configuration space measurement while Tr/k for Fourier space 
measurement. The diamonds/filled squares have been offset in the 
horizontally negative/positive direction for clarity. 

the growth rate is assumed to be a constant when modi- 
fied gravity theories are tested, the ratio ^(*)/^('') cannot 
be used in a simple way for this purpose. Third, the 
quadrupole-to-monopole ratio in Fourier space P^^/") = 
P2/P0 (equation ([6|)) gives almost the same value of /3 
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analysis. The spatial velocity bias should exist, which is 
coupled with the nonlinear stochastic bias discussed in 
the text. 

On the other hand, it is known that the clustering of 
small dark matter halos depends not only on their mass 
but also o n. their assembly hi story, so called the a ssembly 
bias fe.g-. lGao et~all I2OO50 . iWang et al.l H20Q7f ) showed 
that old and low-mass halos that are preferentially associ- 
ated with a high density field are more strongly clustered 
than young halos with the same mass, and consequently 
have higher velocities. Be sides, the stocha sticity between 
halos and dark matter is (|Dekel fc Lahavl lQQQ) can be a 
source of the systematic erro rs in the reconstruc tion. 
Using a method introduced b ylTaruva &: Sutol (I2000D and 
applied to simulation data bv lYoshikawa et all ()200l[ ). we 
have found that both the nonlinearity and the stochas- 
ticity of small halos become larger than massive halos. 
Particularly the stochastic bias monotonically increases 
as the mass of hal o s decr eases, as was found in real space 
bv iHamaus et al.l (|2010D using the two-point statistics. 
Thus the strong scale dependence of /3 for low mass ha- 
los could be caused by the assembly and/or nonlinear 
stochastic bias. However fortunately, the scale depen- 
dence of the measured /3 weakens with the increase of 
halo mass. For massive halos with b > 1.5, the measured 
P approaches the constant predicted by Kaiser's formula 
on scales k < 0.08 h Mpc"^ or r > 30 h'^ Mpc. 

Because the analysis of redshift-space distortions is 
powerful to investigate not only properties of dark en- 
ergy but also modified gravity theories, it will keep on 
playing a key role in ongoing and upcoming large redshift 
surveys, such as BOSS, FMOS, HETDEX, WiggleZ, and 
BigBOSS. Galaxies targeted by the BOSS survey are lu- 
minous red galaxies, which reside in massive halos. In 
this work we demonstrate the /? value can be well re- 
constructed with a redshift distortion analysis of LRGs. 
On the other hand, one of the samples targeted by the 
BigBOSS, for example, is that of emission-line galaxies, 
which reside in halos with a broad range of halo mass. 
One needs to be careful in using such a sample to con- 
strain the growth rate from the redshift distortion, be- 
cause it can be a scale-dependent function. While one 
might be able to obtain a result consistent (or inconsis- 
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tent) with the prediction from general relativity, it could 
be just a coincidence after the scale-dependent growth 
rate is averaged over some separation or wavenumber 
ranges. We will use semi-analytical modeling or a halo 
occupation model to investigate this issue in future work. 

Recently, an in teresting method was proposed by 
iSeliak et al.l ()2009D to suppress the shot noise in power 
spectrum measurement. They considered an optimal 
weighting function f{M) in measuring the galaxy over- 
density, where they give higher weights on higher mass 
halos. Compared with our results presented here, such 
a mass weighting scheme is useful not only for suppress- 
ing the shot noise but also obtaining the true value of 
the growth rate. This scheme can be naturally incor- 
porated into our method and such a study will be pre- 
sented in our future paper. Another theoretical improve- 
ment to be applied for observation is evading the cosmic 
variance limit, which is one of the most important tasks 
for precise measurement of the redsh ift-space distortions, 
as we have already seen above. iMcDonald fc SeliaM 
(|2009t) showed using multiple tracers of density with dif- 
ferent biases suppresses the noise for measurement of 
/3 on large scales dramatically comp ared to the tradi- 
tional single tracer in ethod (see also IWhite et al.l 120091 : 
iGil-Marm et al.ll2010[ ). But the different scale-dependent 
properties of /3 for different halo masses found in Figure 
4 imply that the real situation might be more complex, 
and realistic models of galaxies must be adopted to in- 
vestigate if the method of multiple tracers works. 
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